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Abstract 

> 

^sO ■ A general theorem on factorization of matrices with polynomial entries is proven 

and it is used to reduce polynomial Darboux matrices to linear ones. Some new 
^vq ' examples of linear Darboux matrices are discussed. 

o ' 
in 

1 Introduction 



In the continuous case, for spectral problems with mxm matrix potentials U polynomial 
in A 

d 



the construction of Darboux transformation 

V^V = Q(x,\)y, ^-Q = UQ-QU ^ = U(x,X)^ (2) 
can be based on the elementary Darboux transformations, for which 

Q(x,X) = (X-a)P(x) + (X-/3)(E-P(x)), P 2 = P. (3) 

Here, E denotes the identity matrix and the projector P is uniquely defined by the null- 
spaces of Q : 

KerQ\x =a = Im P = M, Ker Q\ x=f3 = Ker P = N 

1 E-mail : musso@fis . uniroma3 . it 
2 E-mail:shabat@fis. uniroma3.it 



If v(x) G M and w(x) G N, then the following equations have to be satisfied: 
d d 

—v(x) — U(x,a)v(x), — w(x) — U(x, 3)w(x). (4) 
ax dx 

The zeroes of the determinant define the eigenvalues X = a, (3 introduced by the Darboux 
transformation (3). In particular, for the Zakharov-Shabat spectral problems (see, e.g., 

U = \a + q(x), q(x) = [a,q{x)\, (5) 

where a is a constant matrix with simple eigenvalues, the function ^ satisfies the equation 
(1) with the potential 

U = U+[Q,a]. 

The case 

p = a, P* = P (6) 
where * means hermitian conjugation, corresponds to the self-adjoint potentials 

U*(x,X) = U(X,x), q* = q. 

Main Theorem in Section 2 3 formulates the necessary and sufficient conditions on the 
iVth degree matrix polynomial A(X) to be represented as the product of iV first order 
polynomials (3), (6). 

In application to Darboux transformations this factorization A = Qn • ■ • Qi gives rise 
the sequence of spectral problems (1) 

= Uj{x, X)*,, = Qj(x, A)* j5 (7) 

related with each other by the transformations (2), (3). Observe that the composition of 
Darboux transformations ^j+i = is also a Darboux transformation since (Cf. (2)) 

Bj,x — Uj+iBj — BjUj, B < ^B n ---B 1 =3- B x — U n+1 B — BXJ\. (8) 

For a scalar differential operator A, the elementary Darboux transformation ^ — > ^ 
is defined by the first order differential operator: 

A^ = X^, * = (£)-/)*, = X^ A= (D - f)A(D - f)-\ 

For the operator A of order N, an analog of Main Theorem yields the relation 

A* = <^>- =( D -f N )...(D- f 2 ){D - /i)*, (9) 
< <Pi, ■ ■ ■ ,<Pn > 

where ipj G Ker A and 

<^i,...,^>= f det( J D fe - 1 (^)), j,k = l,...,l, D=± 



3 An analog of this theorem was recently used by [7] in order to construct an approximate solution of 
Riemann-Hilbert problem. See also [6]. 
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The proof of (9) is obtained by induction (Cf. Main Theorem) based upon the identity- 
valid for any smooth functions ifij : 

<(Pl,...,(Pm>=<Pl<<f>2,---,<Pm>, <Pj = ( D - fl)<Pj, fl=Dlog(p 1 . 

In eq. (9) the factors D — fj are irreducible (they are first order differential operators). 
Such property appears to be crucial in the theory of dressing chains [2], [3], [4], where 
the functions fj are used as basic dynamical variables. On the other hand, the Darboux 
transformation (3) can be factorized into the composition of two more elementary ones 
with null spaces Ker P and Im P, respectively. This factorization reformulated in terms 
of the corresponding matrices B 1 and B 2 gives rise to the following equations: 

\I/(m + 1, n) — Bi{m, n)^(m, n), \I/(m, n + 1) = B 2 (m, n)ty(m, n), 

(m + 1, n + 1) = Q(m, n)^(m, n), Q = B 1 o B 2 , 

where by definition 

Q = B x oB 2 Q(m, n) — 5i(m, n + l)B 2 (m, n) — B 2 (m+ l,n)B 1 (m, n). (10) 

Equations (10) show that, unlike the continuous case (8), the Darboux transformations 
and the original discrete spectral problem have the same nature and it is possible to 
replace one with the other. 

Thus, the problem of factorization of (3), considered in Section 3, brings up the discrete 
spectral problems with potentials linear in the spectral parameter A. In the examples, 
we will restrict ourself to the simplest 2x2 case but without the assumption that Q is 
factorized through orthogonal projectors (see (3)). In our considerations, the continuous 
independent variable x will play the role of a parameter . 

In the last subsection we consider solutions of the Yang-Baxter equation linear in the 
spectral parameter A. 

2 Main Theorem 

In order to prove that the Darboux transformations (3) are generic ones, it suflicies to 
prove that the product of jV matrices (3) yields an iV-th degree polynomial in A of a 
general form. We will use the characteristic property of the polynomials (3) as follows 

Q a p d ^ f (A - a)P + (A -/3)(E- P) => Q aP Q Pa = (A - a)(X - (3)E. (11) 
In case (6), we have 

Q(A)*Q(A) = (\-a){X-a)E 

, and such property holds, obviously, for the product of matrices (3) with the orthogonal 
projectors P* = P. 
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Theorem. Let A(X) be an m x m matrix with polynomial entries in A satisfying the 
conditions: 

N 

A(\yA(\) = H(\-a j )(\-a j )E, (12) 

lim^oo^ = E. (13) 
Then A(X) can be factorized as 

N 

aw = n (( a - + ( a - %)(^ - p ^)) . (14) 

J'=l 

where the P aj . are Hermitian projectors. 

Proof: We proceed by induction on N. UN — 0, then from (13) and the Liouville 
theorem it follows that A(X) = E. 

Now we show that if the statement of our theorem holds for N — 1, it holds for TV as 
well. From (12) and (13) it follows: 

N 

det(A(X)) = n(A - a,-)* (A - . (15) 

j'=i 



We select a zero of det(A(A)); from (15) it follows: 

(16) 



dim(Ker(A(<x,))) = A^-, 



dim(Ker(A(<x,))) = m — k 



,r 



Let | Vi)^=i and ju^)^^ be an orthonormal basis for Ker (A(otj)) and Ker(A(<x,)), respec- 
tively. We denote by P aj and P^ the Hermitian projectors 

^ = EN)<4 (i7) 
i=i 

P a . = ]T \wi){wi\. (18) 
i=i 

Let us prove that 

P Sj = (E-P aj ). (19) 
To prove (19) is equivalent to proving that C m admits an orthogonal decomposition 

C m = Kev(A( aj )) © Ker(A(aj)). (20) 

From (12) it follows: 

A(a,)A(^T = 0, (21) 

that is 

Im(A(^)*) C Ker(A( aj )); (22) 
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but from (16) we have: 



whence it follows: 



dim(Im(A(a!j)*)) = m — dim(Ker(A(<x,))) = kj, 
dim(Ker(A(<x,))) = kj, 

lm(A(a i y)=Ker(A(a 3 )). 



Equation (23) is equivalent to 

Ker(A(^)) = lm(A( aj y) 

whence it follows 

C m = Ker(A( aj )) © lm(A( aj )*) = Ker(A(aj)) © Ker(A(a j ))- 
Now the matrix 

A(X) [(A - aj)P a . + (A - aj )(E - P a .j 
vanishes at A = ctj and A = otj. Consequently, we have: 

A(\) [(A - aj)P a . + (A - aj )(E - P a .)] = (A - a,) (A - a,)i(A) 

for a certain matrix A(X). Then 

A(X) = A(\) [(A - aj )P a . + (A - otj){E - P a .j 

Since P a is Hermitian, it follows that A(X) satisfies 

N 



Moreover, 



A(\yA(\) = l[(\-a k )(\-a k )E. 

k+3 



lim A -^oo4^ = lim A ^ooT^r-r = E. 



In the case of distinct zeroes otj ^ a k , the orthogonal projectors P a . in ('. 
uniquely defined by the subspaces 

M J d ^ f KerA(A)| A=Qj 
and an ordering of zeroes a>\, . . . , a at. Indeed, for instance, for = 2, we have 

A(X) = Qi(A)Q 2 (A), Qj(X) = (A - aj )P aj + (A - a 3 )(E - P a .), j = 1, 2 
with P* = P a and 

ImP Q2 = Ker(£ - P a2 ) = KerQ 2 (a 2 ) = Ker Q 1 {a 2 )Q 2 {a 2 ) = M 2 



5 



, since det Qi(a 2 ) ^ 0. Similarly 

ImP Ql =Q 2 (a 1 )M 1 , since M 1 = Ker(E - P aj )Q 2 (ai). (31) 

In order to construct a dressing chain, we have to factorize (see the next section) the 
basic Darboux transformation (11) further: 

Qap(X) = (A - a)P a(3 + (A - (3)(E - P af3 ) = B 2 (X)B 1 (X) (32) 

where these two factors Bj(X) = crjX + bj are linear in A and 

det Q a/3 (X) = (A - a) k {X - p) m ~ k det B x = (A - a) k , det B 2 = (A - /3)" 1 "*. 

Evidently, 

Ker = Im S 2 (a) = KerP a/3 , (33) 

since, due to (32) and the fact that det Bi(a) ^ 0, we have 

Q a p(a) = {a- 0)(E - P a p) = B 2 (a)Bi(a) Im^a) = lm(E - P a/3 ) = KerP a/3 . 

3 2x2 case 

In this section we consider some interesting, from our point of view, classes of 2 x 2 matrices 
linear in A related with (3). First, let us consider the chain of matrices Bj, j G Z : 

Bj = ^ \^ 2 ^ ^ , hj = ^Uj Vj+1 - (3j for j odd, (34) 
Bj = ^ \ 2 +hj ) ' h i = \ v i u J+i ~ & for 3 even - ( 35 ) 

Thus 

det Bj = A — f3 j for all j 
and (Cf (33)) it easy to see that 

Ker Bi | A=/3l = Im B 2 \ x=f 3 2 , Ker B 2 \ x =p 2 = Im B 3 | A=A 

Moreover, these matrices satisfy the consistency conditions (10) (see Lemma below) and 

B 2 B X = X - ft + f3 12 P 12 , B 3 B 2 = X - f3 2 + (3 23 P 23 , /% d =i f ft - ft, 
where Pfj = Pij are projectors. In particular, 

o r> ( \v 2 u X3 \u 13 \ def 

,n 73 / (^ 2 3 + jU 2 v 42 ^u 3 (p 32 + \u 2 v 24 ) \ def 

P23-T23 — 1 , 1 , „, , % 
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^42 4^3^24 



Thus these matrices provide us with a factorization of elementary Darboux transforma- 
tions (3). Observe that, for odd i, the product BjB i: where j — i + 1, depends only on 
for z odd and on the difference v i+2 — v i+ i and u i+1 for i even. 
In comparison with the factorization defined by main Theorem, the entries of matrices 
(34) and (35) are expressed directly in terms of the potentials 

U '=(v } -x) =Xo + * 

of the 2x2 Zakharov-Shabat spectral problem 

$\ - A^ 1 = u^ 2 , Vl + \V 2 = vV 1 , or V X = UV. (36) 

Thus, the substitution of matrices (34) and (35) into the equations (8) yields the dressing 
chain for the spectral problem (36): 

(D + 2/3j)v j+1 + 2vj = ^UjV 2 +1 , (D - 2^)uj - 2u j+1 = -^u 2 jVj+1 , (37) 

(D - 2/3j)u j+1 - 2u 3 = -^u 2 +lVj , (D + 2fa) Vj + 2v j+1 = ^u j+1 v 2 . (38) 
The first case corresponds to (34) and the second one to (35). 

Remark. In the case of the general Zakharov-Shabat spectral problem (5) one could 
consider the following Darboux transformations (8) 

Bj = Ao"j + bj, Uj = Act + qy 

Collecting quadratic and linear terms in A in (8) one finds 

[a, Gj] = 0, a j:X = [a, Aj] + q j+1 aj - 

That implies that o^ x = since the potentials qj are off-diagonal in the basis where the 
matrix a is diagonal. Thus, the matrices Uj should be x— independent and could be any 
diagonal matrix in that basis. 

Coming back to the spectral problem (36) we will show how, starting with (ui,vi), 
one construct a chain of Darboux matrices Bj, transforming the equation (36) with the 
potential Uj into one with the potential Uj+i, using the solution to the corresponding 
spectral problem Riccati equation: 

p x = 2\p - Vjp 2 +Uj, p d = — . (39) 
This should highlight a similarity with the dressing chain considered in [2] . 
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Proposition. There is a one-to-one correspondence between the matrices (34), (35) and 
solutions of the Riccati equation (39) with A = (3j. 

Proof: At the first step we have to find the potentials u 2 and v 2 . The function v 2 is 
defined by (37) as solution of ODE 

(D + 2/3i)u2 + 2ui = l -u x v\ 

which is reduced by rescaling v 2 = 2p to (39) with A = f3\. For the function u 2 , we obtain 
now the exact formula ^ 

2u 2 = u 1)X - 2(3x111 + -ju\v 2 . 

At the next step, we find w 3 = 2p _1 by solving (39) with A = /3 2 - It satisfies the 
corresponding equation (38): 

(D - 2/3 2 )u 3 - 2u 2 = -u\v 2 

while v 3 is defined by the second equation in (38). Thus in both cases the system of 
equations (37) and (38) are reduced to the Riccati equation (39). I 

Consistency problem. We will ignore now the dependence on the continuous indepen- 
dent variable x and consider more general (as comparised with (34), (35)) matrices linear 
in A and satisfying conditions analogous to (10). Having applications in mind, we will 
denote one of the matrices in (10) as A(n, A) and consider the second one as a Darboux 
transformation for the discrete spectral problem defined by A : 

y n+1 = A(n,\)V n , A=^ X + ° °y * = # n+1 = i(n,A)*. (40) 

If the spectral problem in (40) is a "good" one, 4 then there should exist nontrivial 
Darboux transformations with 

*-(»*-*) (41) 

and the consistency equations T 1 oT 2 = T 2 o T±, where = AV, T 2 ^ = B^> (Cf. (10)) 
yield 

(A + c)e + ga = e'(\ + c) + fb, (A + c)f + a(h - A) = e'a + fa 
eb + ag = (A + c)g' + b(ti - A), a(h - X)e + fb = a(ti - A) + g'a 

where B' d = T^B) and A d =T 2 (A). We see that 

g' = T 1 (g) = b, a = T 2 (a) = f\ e' = e, a = a 

and thus the consistency equations are reduced to 

| b(h' + c)=et> + ag, e(c - c) = f'g[ - fg ^ 
[ a(h + c) = ea + otf, a(h' — h) — ab — ab 

It is easy to verify now 

4 The particular case of A with a = is related with Toda chain. 
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Lemma. The matrices Bj and Bj +1 defined by (34), (35) satisfy the consistency equa- 
tions (42) for any j. 

Generally we find that (Cf. (10)) 

( e(\-5 + h!) ab> \ 

AoB ~ \ btf a(tf-\ + P) ) 

where 

fit = c+h = c+ti, ab = a(c+(3), fg = e(h-5), det A = a(A-/?i), det B = e((3 2 -X) 

and thus our Lemma corresponds to the case e + a = 0. 

The Darboux matrix B in (40) is reciprocal to A and (unlike the continuous case) one 
can replace it with Ao B. This property allows us to use the matrix A as a starting point 
in the spectral problem (40) with a = [3]. 

An i?-matrix interlude. While working out examples for the above factorization the- 
orem, we noticed that some of them were also solutions of the quantum Yang-Baxter 
equation 

#12 (A - /i)i?i 3 (A) J R 23 (/i) = #23(^)#i3(A) J Ri 2 (A - fi). (43) 

Here R 12 = R<£> Id, R 23 = Id® R and R 13 = (Id <g> P)R 12 (Id <g> P) acts trivially in the 
second of the three spaces, P is the usual permutation operator: 

P(u ® v) ^ v ® u. 

Thus in the simplest 2x2 case, (43) is an equation for a 4 x 4 matrix R and a well known 
solution is given by 



R(X) 



( X + ir] ^ 

X ir] 

27/ A 

V X + ir] J 



l[(X + i v )P + (X-ir,)(l-P)}. (44) 



and corresponds to the XXX— spin chain. Since P is a projector, we can express it in 
the form: 

P = MiM, wh e r e 2<a| = <6j = (0, 1, -1, 0) (45) 
(a\b) 

We found that there exist two other possible choices of (a\ and (b\ in (45) such that an 
i?-matrix of the form 

R(\) = (\ + i7j)P + (\-irj)(l-P) (46) 
satisfies the Yang-Baxter equation, namely: 

(a\ = (a~\ 1, 1, a), (b\ = (/T\ 1, 1, (3), 

(a\ = (—1, a, (3, —13a), (b\ = ((3a, a, (3, 1), 
with a, (3 being two real parameters. 
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Having two new solutions of the Yang-Baxter equation, we can construct the associ- 
ated Lax matrices and, finally, two monodromy matrices out of them. Then, taking the 
traces of the monodromy matrices we obtain two sets of commuting quantum observables. 
Unfortunately the elements of each set are functionally dependent. Hence, these two new 
solutions do not give rise to new quantum integrable systems. 

Further investigations on Quantum Yang-Baxter equation showed that, if we gen- 
eralize the i?-matrix form (46) dropping the normalization condition (13), then other 
solutions, linear in the spectral parameter A, can be found. We notice that R can then 
be expressed in the form R = XA + B, where A and B must be solutions to the two- 
dimensional constant quantum Yang-Baxter equation; all such solutions have been found 
by J. Hietarinta [5]. 

An example of solution that we believe could be of some interest is the following one: 

A \ / — 77 A \ 

! _ J_ A - r, 

0' " X 2 -rf A — 77 
77 / VAOO-77/ 

Work is in progress on this topic. 

Acknowledgements We would like to thank Decio Levi and Antonio Degasperis for 
useful discussions and the interest shown in this paper. 

This work was partially supported by Russian Foundation for Basic Research, grants 
No. 04-01-00403 and Nsh 1716.2003.1. 

References 

[1] V.Zakharov and A.Shabat, "A Scheme for integration the nonlinear equations of 
mathematical physics by the method of the inverse scattering problem, II", Fund. 
Anal. Appl. 13: 166-74, 1979. 

[2] A. P. Veselov and A.B. Shabat "Dressing chain and spectral theory of Schrodinger 
operator", Funkts. Anal. Pnlozen., 27(2): 1-21, 1993. 

[3] A.B. Shabat, "Third version of the dressing method" Theor. Math. Phys., 121(1): 
1397-408, 1999. 

[4] V.E.Adler, V.G. Marikhin, A.B. Shabat, "Canonical Backlund transformations and 
Lagrangian chains," Theor. Math. Phys., 129(2): 163-183, 2001. 

[5] J. Hietarinta, "Solving the two-dimensional constant quantum Yang-Baxter equa- 
tion" J. Math. Phys., 34(5): 1725-1756, 1993. 

[6] I.T.Habibullin, A.Shagalov, "Numerical realization of the IS," Theor. Math. Ph., 
83(3): 565-73 (323-33), 1990. 

[7] A.V.Domrin, "On a version of the Riemann problem in the theory of integrable 
systems," Dokl. Ross. Akad. Nauk, 390: 301-303, 2003. 



R{\) 



( r] 

X + n 

A + 77 

V A 



10 



